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1. Introduction 

GL-QCD is a QCD-like theory in which the gauge group SU (3) of strong interactions is replaced by 
the exceptional Lie group G 2 . The theory is a gauge theory with fermionic baryons and fundamental 
quarks [1,2] and it can be simulated without sign problem at finite density and temperature. Unlike 
other QCD-like theories such as adjoint QCD or two color QCD, for example, its properties in 
the quenched case are very si mi lar to those of QCD [3-8]. Although the center of G 2 is trivial, it 
shows a first order deconfinement transition which has quite interesting implications for the role of 
the center symmetry in QCD as reviewed in [9] for Lattice 2012. 

In this contribution we briefly summarize our previous results [10] for the hadron spectrum in 
the vacuum and the phase structure at zero temperature as seen in the quark density from lattice 
simulations on rather small 8 3 x 16 lattices in Sections 2 and 3. The observed structures in the 
density over the chemical potential can thereby be related to a corresponding hierarchy of mass 
scales in the baryon spectrum. In particular, one observes thresholds in the baryon density at 
values of the chemical potential that correspond to the pseudo-Goldstone scalar diquark scale, 
an intermediate pseudo-scalar and vector diquark scale, and roughly the fermionic baryon scale 
set by the GL-nucIcons and delta baryons. In Section 4 we furthermore present results from our 
latest simulations in which we have accumulated evidence for a zero-temperature first-order phase 
transition to what might be called G2-nuclear matter. 

The most important difference to QCD is the existence of diquarks in the hadronic spectrum. 
The lightest two diquark states arc the pseudo-Goldstone bosons of chiral symmetry breaking, and 
their quantum numbers differ only by charge conjugation. In order to investigate chiral symmetry 
breaking in detail we add diquark source terms to disentangle the charge conjugation symmetry. 
Therefore we have to use Majorana fermions in the simulations. In Section 5 we review the formu¬ 
lation of G2-QCD in terms of Majorana fermions and discuss the chiral properties in the presence 
of diquark sources in more detail. First preliminary results of our simulations at finite temperature 
and density arc shown in Sections 6 and 7. 

2. Chiral symmetry and baryon number in G 2 -QCD 

The Euclidean action of Nf = 1 flavour LL-QCD with quark chemical potential /.( for baryons reads 

S = J d 4 .x | - ^ tr Up v F ‘ uv + ¥ D [A , m, fit} ¥ 

D\A,m,n] =y^(d^-gA^-m + y 0 fi, 

where the gauge group is the exceptional Lie group G 2 , and y u = y u arc the Euclidean /-matrices. 
The fundamental representations of G 2 arc 7-dimensional and 14-dimensional, the latter coinciding 
with the adjoint representation. Since G 2 is a subgroup of 50(7), all representations arc real. The 
Dirac operator satisfies 

000+75 = 75 £>(-AU) and D(fi)*T = TD(n*) (2.2) 

with T = C75, T* T = —1, T 1 = T 1 and charge conjugation matrix C. If such a unitary opera¬ 
tor T exists then the eigenvalues of the Dirac operator come in complex conjugate pairs, all real 


with 
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Figure 1: Pattern of chiral symmetry breaking in G 2 -QCD. 


eigenvalues are doubly degenerate [11,12] and thus 

detD[A,m,fi\ > 0 for gGl. (2.3) 

GL-QCD with a single Dirac flavour possesses an extended chiral symmetry [11] compared to 
QCD. The action is invariant under the SO(2)\j vector transformations and the usual axial trans¬ 
formations leading to a U (2) symmetry group, in agreement with the results in [1]. Following the 
same arguments as in QCD it is expected that the axial (7(1) is broken by the axial anomaly such 
that only a 517(2) x Z(2 )b chiral symmetry remains. In the presence of a non-vanishing Dirac 
mass term (or a non-vanishing chiral condensate) the theory is no longer invariant under the axial 
transformations. Therefore the non-anomalous chiral symmetry is expected to be broken explicitly 
(or spontaneously) to its maximal vector subgroup, 

517(2) ®Z(2)b^50(2)v®Z( 2) b , (2.4) 

The remaining chiral symmetry at finite baryon chemical potential is the same as in QCD, 

57/(2) <g>Z(2) B ^7/(l) B (2.5) 

The final pattern of chiral symmetry breaking of G 2 -QCD is shown in Figure [F]. If chiral symmetry 
is spontaneously broken, the axial chiral multiplet becomes massless, according to the Goldstone 
theorem. The following operators generate the two Goldstone bosons: 

d(0 ++ ) = 'f /C y 5 'f' - x f'y 5 T /C and r/(0 + ") ='f /C 7 5 x f / + x f / 7 5 'f /C . (2.6) 

They have quark number n q = 2 and hence also carry a baryon number of /?b = 2/3, if baryon 
number counts the difference of quarks and anti-quarks per f/Vnuclcon as in QCD. The Goldstone 
bosons in CL-QCD are scalar diquarks instead of pseudoscalar mesons as in QCD. As long as we 
do not introduce diquark sources, r/(0 ++ ) and r/(0 1 ) have the same mass. In Section 5 below we 
also introduce the corresponding diquark source terms, however, in order to disentangle states with 
opposite charge quantum number and investigate chiral symmetry breaking more closely. 
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3. Spectroscopy and the phase diagram at zero temperature 

The possible quark and gluon content of (colorless) bound states is determined by the tensor prod¬ 
ucts of the appropiate representations of the gauge group Gi. Quarks in Gi transform under the 
7-dimensional fundamental representation, gluons under the 14-dimensional fundamental (and at 
the same time adjoint) representation. An overview over possible bound states can be found in [10]. 
We expect to find bound states for every integer quark number n q . Mesons have n q = 0, diquarks 
n q = 2, and nucleons n q = 3. In addition, there are more exotic bound states of gluons and quarks, 
for example a hybrid with n q = 1. In the following we give an overview over the bound states 
considered here, where u and d denote flavours of Dirac fermions. For the nucleons N and the pion 
n we make use of the partially quenched approximation in our one-flavour simulations. Table |] 
shows bound states that are also present in QCD while Table [2] shows the diquarks. In all tables 


Table 1: Bound states of G 2 -QCD with 2 flavours for baryon number «b = 0 (left) and baryon number 
= 1, he. quark number n q = 3 (right). For details see text. 


Name 

0 

T 

J 

p 

c 

71 

iiy 5 d 

SASS 

0 

- 

+ 

V 

iiy 5 u 

SASS 

0 

- 

+ 


Name 

0 

T 

J 

P 

C 

N 

T abc (u^y 5 d b )u c 

SAAA 

1/2 

± 

± 

A 

T abc (u%Y fl u b )u c 

SSAS 

3/2 

± 

± 


Name 

0 

T 

J 

p 

c 

r/(0 ++ ) 

u c Y5U + c.c. 

SASS 

0 

+ 

+ 

d( 0 +-) 

ii C 75 u — c.c. 

SASS 

0 

+ 

- 

d( 0 -+) 

uru + c.c. 

SASS 

0 

- 

+ 

o_ 

“e 

U C ll — c.c. 

SASS 

0 

- 

- 

d(l++) 

vry^d ' — d c Y b u + c.c. 

SSSA 

1 

+ 

+ 

d( 1 + ") 

u C Ynd — d c Yn 11 ~ c - c ■ 

SSSA 

1 

+ 

- 

d( 1 - + ) 

ti C Y5Ytid - d c Y5Yti u + c.c. 

SSSA 

1 

- 

+ 

d( I") 

d C Y7 Yu d ~ d c Y5Yn u ~ c.c. 

SSSA 

1 

- 

- 


Table 2: Bound states with baryon number «b = 2/3, i.e. quark number n q = 2. 


0 is the interpolating operator used to extract the mass in simulations, the string T represents the 
behaviour of the wave function under change of position, spin, colour and flavour (S stands for 
symmetric, A for anti-symmetric), and J, P, C are the spin, parity and charge conjugation quantum 
numbers. The difference between the rj and the diquark correlation function is only the discon¬ 
nected contribution. Therefore, the diquark with positive parity has the same mass a the pion with 
negative parity, m d ( 0 +) = m,(o~)- I n [ 10 ] it i s shown that for every diquark there is a flavour non¬ 
singlet meson with the same mass but opposite parity. In the following we discuss two different 
ensembles with lattice parameters as listed in Table |3| A physical scale is set by the proton mass, 
= 938 MeV. The mass spectrum for both ensembles is shown in Figure |2| In the heavy en¬ 
semble the diquark masses and all parity even and odd states arc almost degenerate. In the light 
ensemble the diquark masses are no longer degenerate. We observe a significant mass splitting 
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Ensemble 

P 

K 

m d( 0+) a 

m^a 

m d( 0 +) [MeV] 

a [fm] 

a -1 [MeV] 

MC 

Heavy 

1.05 

0.147 

0.59(2) 

1.70(9) 

326 

0.357(33) 

552(50) 

7K 

Light 

0.96 

0.159 

0.43(2) 

1.63(13) 

247 

0.343(45) 

575(75) 

5K 


Table 3: Parameters for the two different ensembles. All results are from 8 3 x 16 lattices. 


between parity even and odd states as well as between scalar and vector diquarks. Especially the 
Goldstone boson becomes the lightest state, with the r\ also being somewhat heavier. For the nu¬ 
cleons we also observe different masses for parity even and odd states and the spin 1/2 and spin 
3/2 representations. In particuclar, we find three clearly different scales in the light spectrum: a 
pseudo-Goldstone scale, an intermediate boson scale set by the remaining diquarks, and the nu¬ 
cleon scale set by the N and A masses. This mass hierarchy of the spectrum seems to be reflected 
in various structures of the quark density at zero temperature which one might thus attribute to 
different bosonic and fermionic phases at finite density, see Figure [ij. With increasing chemical 
potential, the quark number density first remains consistent with zero until it very quickly rises 
to a very small but nonzero value. When we compare the critical chemical potential fd c for this 
onset transition to the mass of the lightest baryon m^o+b the pseudo-Goldstone 0 + diquark in our 
case, we find that numerically very good agreement with the expectation from the Silver Blaze 
property, i.e. = m d ( 0 +\/2. The ground state changes from the vacuum to a finite-density ground 
state only when the quark chemical potential reaches the mass of the lightest baryon divided by its 
quark number so that the corresponding excitation energy vanishes. For bosonic excitations one 
might expect Bose-Einstein-condensation in a continuous second-order quantum phase transition 
at jU c = m d({) \ j/2, without binding energy, and our data is certainly consistent with that. For larger 
values of the chemical potential plateaus develop where the quark number density remains almost 
constant. Especially in the light ensemble, the step towards the second plateau conicides with the 
mass of the heavier bosonic diquark states divided by their quark number. It appears that the two 
bosonic baryon mass scales are not sufficiently separated from each other to resolve these two 
distinct transitions in the heavier ensemble. 

At around afi = 0.6 for the heavy ensemble and afl = 0.55 for the light ensemble the quark 



Figure 2: Mass spectrum for the heavy (left) and light (right) ensemble in G 2 -QCD. 
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Figure 3: Quark number density for the heavy (left) and light (right) ensemble in GVQCD at finite density 
and zero temperature. 


number density starts increasing again and no further plateau is observed. This transition appears 
to coincide with the mass scale of nucleon and A divided by three. In both ensembles the general 
pattern thus seems to be that the various transitions in the quark number density are related to 
the various baryon masses in units of their quark number. While for bosonic baryons the density 
bends towards the zero-axis with plateaus forming after each transition, at the scale of the fermionic 
baryons the quark number density is convex, as seen most clearly in the light ensemble with better 
separation of scales, and continues to further increase with increasing chemical potential until 
saturation sets in, eventually, when the lattice starts to get filled with the maximum number of 
quarks per site, i.e. at el'll q = 14 here. This clearly is a lattice ai tifact beyond the range of any 
hadronic interpretation of the density and it is therefore not shown here again, see [10]. 


/i in MeV 

180 200 220 240 260 



afl 


Figure 4: Quark number density for the light ensemble in the vicinity of the transition. 
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4. Evidence of a first order nuclear matter transition 

In both ensembles we observe a strong transition at a/i ~ 0.52 (heavy ensemble) and a/i ~ 0.38 
(light ensemble) that does not appear to cone pond to any of our spectroscopic states. In Figure |] 
we show the vicinity of this transition in the light ensemble in more detail. The quark number 
density rises between afi = 0.36 and afi = 0.40 from a lower value n q « 0.010 to a higher value 
n q ~ 0.025. In Figure || we show the quark number density as a function of Monte-Carlo time and 
observe tunneling between these two states. This might indicate that there is a first-order phase 
transition at afi ~ 0.38 in the phase diagram at zero temperature. Whether this phase transition is 
indeed the analogue of the liquid-gas transition of nuclear matter as expected in QCD remains to 
be shown by further simulations. If this is the case, then either the binding energy per nucleon is 
comparatively large or the masses of nucleon and A change with density in the regime of the finite 
bosonic baryon density in the ground state before this transition which is not possible in QCD. 




0 500 1000 1500 2000 2500 3000 3500 

Monte Carlo time 


0 500 1000 1500 2000 2500 3000 3500 

Monte Carlo time 


0 500 1000 1500 2000 2500 3000 3500 
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Figure 5: Quark number density as a function of Monte-Carlo time for different values of chemical potential. 
From left to right in the upper row afi = 0.36, a Li = 0.37 and all = 0.38 and in the lower row a/i = 0.39, 
afi = 0.40 and a/i = 0.41. 


5. Simulations of G2-QCD with Majorana fermions 

In the present section we discuss the introduction of diquark sources in G 2 -QCD for two reasons: 
First, on larger lattices and especially for values of the chemical potential in the vicinity of the 
first order transition the simulations become more and more expensive. An obvious reason for this 
might be the presence of very light diquark excitations in the simulations in this region of the phase 
diagram. With the introduction of diquark source terms simulations should become more feasible. 
Secondly, we would eventually like to resolve the complete diquark spectrum in order to investigate 
chiral symmetry breaking at finite density. Any diquark source necessarily consists of an operator 
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with quark number « q = 2 and therefore contains a pair of charge conjugated Dirac spinors 'F for 
the quarks. In order to integrate over the fermion fields in the path integral, however, one needs 
a bilinear expression in V P and V P. Unlike two-color QCD, however, for a single flavor of Dirac 
fermions it is then necessary to introduce corresponding Majorana fermions in the simulations. 
Since the gauge fields in G 2 satisfy Aj t = -Aa = —A " T a it is possible to rewrite the matter part of 
the action in (^1]) as a sum over 2 for /./ = 0 at first uncoupled Majorana spinors A = fix , tj), 

S^,A]=Jd*x'¥(r(d tl -gAi 1 )-m)'¥ = Jd A xl gA^fi)-m) A (5.1) 

Here A obeys the Majorana condition A c = CA T = A, A c = —A T C 1 = A, and it is related to the 
Dirac spinor as'P = ^-|-ir],'i / = ^ — ifj, X F C = X ~ i V and 'F C = X + i V • The (baryon) chemical 
potential /./ is an off-diagonal term in Majorana flavour space such that 

2 z? = x ¥D(/n,ju) x ¥ = AM(m,ju) A with 

D(m,jtt) =0 - m + y (l /./ and = (0 - m)<7 0 - /i/oGb 


where Oq = 1 and the Pauli matrix 02 act on the 2 Majorana flavours X and tj which are now 
coupled to each other by the offdiagonal /./ovtcrm. Integration in the path integral over the Majo¬ 
rana fermions leads to the Pfaffian instead of the fermion determinant and we can proove that the 
Pfaffian is positive as expected. 

It is possible to introduce two different diquark sources in G 2 -QCD, one for the scalar ( 2 z? rs ) 
and one for the pseudoscalar diquarks ( 2 z?i), 

^y 5 (J) = \ , ^i(J) = \ + . (5.3) 

In the Majorana decomposition these terms read (J = J\ + iU- J — J\ + 1 / 2 ) 

=i (fl) (a -a) rs (i?) = iA + ' ,2<rj) ' 

= (fl) (-A -a) (j,) = A {J '* - hc ' )x ■ 

The Lagrange density for the matter part of the theory is then given by 


££ = \ [00 —m-m 5 y 5 )c7o —P7otf 2 + (i/t75 —T 2 )cti + (/i+ 1 / 275 ) 03 ] A. (5.5) 


Similar to the case without diquark sources one can show that the Pfaffian is real if J 2 = di = 0, but it 
is not necessarily positive any longer. Nevertheless, we expect that for small values of the chemical 
potential the sign problem is not present and this expectation is confirmed in our simulations. The 
first derivatives of the partition function with respect to J, J and m define the chiral and diquark 
condensates, 


£ = 
£1 = 

I 5 = 


1 d\n{Z(m,J,J)) 

V dm 

1 dln(Z(m,JJ)) 

V dJ 

1 dln(Z(m,J,J)) 


dJ 


(t a t) = (xx + fm), 

(y Cx ¥ + c.c.^ = (AA-TjTj), 

^'P c y 5 'P + c.c.^ =i(A75fl + fl75A), 


(5.6) 
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Operator 

Parameter 

Ga.i 

Oa.3 

Ov ,2 

Goldstone bosons 

Massive state 

Aa 

m 

X 

X 


<i( 0 ++ ) ,r/( 0 + ) 

/( 0 ++ ) 

A 75 0 i A 


X 


X 

d( 0 ++), /( 0 ++) 

d( 0 +-) 

A7503A 

k 


X 

X 

d( 0 +-),/( 0 ++ ) 

</( 0 ++ ) 

A75A 

m 5 

X 

X 


</( 0 “ + ), rf( 0 “) 

T ?(o~ + ) 

A01A 

Ji 

X 


X 

d(0~ + ) , 77 (0 h ) 

<7(0 ) 

A 03 A 

Jl 


X 

X 

d( 0-),T7(0-+) 

d( 0 “ + ) 

A 7 o 0 2 A 

M 

X 

X 


- 

- 


Table 4: The table shows the transformation behaviour of bilinears under the chiral transformations, / 
means invariant, x not invariant. In the last two columns the corresponding Goldstone bosons and the 
massive states are shown. 

that we investigate in the following. For a single Dirac flavour the chiral symmetry is SU ( 2 ) l = r *. 
The generators for the symmetry transformations are given by Pauli matrices 

r v = l®a 2 , r A = <8 >{cti,<t 3 }, (5.7) 

and the chiral transformations read 

0 A ,i A =e iays ffl A and A -> A e iars 

0 A 3 A= e iarsC73 A and A -> (5.8) 

0 V2 A =e xac -l and X^X e -' xaa -. 

Possible bilinear bound states for a single Dirac flavour are 

d(0 + ~) =Xy 5 ail = xm , d( 0 ++ ) = A 75 <73A = XYsX ~ VYsV , 

d(0—) =1 <t 1 X = XV, d(0~ + ) = Xa 3 X=xx-VV, (5-9) 

/(0++) =XX=XX + VTU ??((r + ) = Ay 5 A = XYsZ + VYs V- 

Table [4] shows their behaviour under the chiral transformations together with the corresponding 
Goldstone bosons. Under the chiral SU (2) these biliners decompose as 

202 = 301, (5.10) 

and we can identify a positive and negative parity triplet, 

(/( 0 ++ ), d( 0 ++), d( 0 +-)) and (17(0 h ), d(0~+), d( 0 “)) . 

Since the negative parity multiplet obtains a contribution to its mass from the chiral 
will set the corresponding sources Ji = J 2 = tn$ = 0. Under a general (infinitesimal) 
formation 5A = i ( 0:7501 + /3 750 2 + 702 ) A the Lagrangian density transforms as 

= 2il ( 750 !(-am + 7 / 2 ) + 7 5 0 3 (-j 8 m- yj x ) + i(aJi +/3/ 2 ) + ift 707 >(/ 30 i - OC 0 2 )) A, 

(5.12) 


(5.11) 

anomaly we 
chiral trans- 
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m 

ha 

Solution 

Symmetry 

Generator 





SU (2) 

T = {75<7l,75^3,a 2 } 


X 


a = (3 = 0 

U( 1) 

T = c 2 



X 

y = 0, aJi = -(3J 2 

U( 1) 

T = 7s(/ 2 CJ| -/i<t 3 ) 


X 

X 

am = yJ 2 ,(3m = —7/1 

U( 1) 

T = 75(/ 2 C7 i - h (7 3 ) + rna 2 

X 



a = (3 = 0 

U{ 1) 

II 

X 

X 


a = (3 = 0 

U( 1) 

II 

X 


X 

a = (3 = 7 = 0 

- 

- 

X 

X 

X 

a = (3 = 7 = 0 

- 

- 


Table 5: Chiral symmetry of Go-QCD in the presence of diquark sources and chemical potential,/ means 
= 0, x means ^ 0. 


and we obtain the following system of equations for the invariance of the Lagrangian 

fia = 0, jU/3=0, am = yf 2 , (3 m = -yh , aJ\ = —(3J 2 . (5.13) 

Possible solutions are shown in Table |5|. The Goldstone bosons are then linear combinations in 
the corresponding multiplet. An illustration of the chiral symmetry breaking is shown in Figure |j. 
At vanishing chemical potential, any linear combination of mass m and diquark sources J\ and J 2 





Figure 6: Chiral symmetry breaking in Go-QCD in the presence of diquark sources. 

breaks the chiral symmetry down to a (/(l) subgroup and baryon number is conserved. The direc¬ 
tion (generator) for the invariant subgroup (7(1) of the SU (2) symmetry is shown in Table [5j Since 
the mass m and the chemical potential ji break the SU (2) in the same direction, baryon number is 
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Figure 7: Chiral condensate, diquark condensate and quark number density as a function of the diquark 
source coupling J\ in the vacuum at T = fl = 0. 


explicilty broken at finite chemical potential by any linear combination of diquark sources. Thus 
in our simulations at zero density we expect no qualitative difference between a fermion mass term 
and a diquark source term. 

6. Results at finite temperature 

Our first aim is to investigate spontaneous chiral symmetry breaking at finite temperature for fi = 0. 
In our simulations with Dirac fermions the corresponding signal in the chiral condensate remained 
rather inconclusive when measured across the deconfinement transition as observed in the Polyakov 
loop [2]. The reason for the weak transition in the chiral condensate was probably the explicit 
breaking due to the Wilson mass. The first simulations with Majorana fermions at zero temper¬ 
ature have been performed on a 8 3 x 16 lattice with /3 = 0.96 and K = 0.151. Compared to the 
ensembles discussed for Dirac fermions before, this corresponds to very heavy quarks (diquarks). 
For vanishing chemical potential the Pfaffian is again positive and simulations can be done. In 
Figure [7] the chiral condensate, the diquark condensate and the quark number density arc shown as 
a function of the diquark source J\. As expected, with increasing diquark source (and fixed mass 
m), the chiral condensate rotates into a diquark condensate. Furthermore, the quark number den¬ 
sity vanishes, indicating an unbroken (7(1) subgroup such that baryon number is still conserved. 
Our simulations at finite temperature have been performed on a 12 3 x 6 lattice with /3 = 0.96 and 
K = 0.156. In Figure ||, the chiral condensate, diquark condensate and the Polyakov loop are shown 
as a function of the inverse gauge coupling /3 for three different values of the diquark source. With 
increasing diquark source the deconfinement transition visible in the Polyakov loop is shifted to 
larger values of /3. This is the same qualitative behaviour as observed when increasing the fermion 
mass. In the chiral condensate and the diquark condensate only a very weak transition is observed 
which becomes more pronounced for smaller diquark masses. For a particular diquark source of 
J\ = 0.20 we have normalized the condensates by their vacuum contribution. The result is shown in 
Figure |9| The transition in the chiral condensate is very weak while it is by a factor of ~ 40 stronger 
in the diquark condensate. We have so far ignored the necessity for additive renormalization in the 
chiral condensate, however. Also the error bars are much smaller for the diquark condensate as 
compared to the chiral condensate. Compared to the situation without diquark sources we observe 
a chiral improvement like for instance in twisted mass QCD where a similar term is added to the 
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Figure 8: Chiral condensate (left), diquark condensate (in units of./ 1 , middle) and Polyakov loop (right) as 
a function of the inverse gauge coupling j3 for three different values of the diquark source /i. 


action in order to improve the chiral properties. The transition temperature agrees with the one 
obtained from the Polyakov loop. A future task in our simulations with Majorana fermions at finite 
temperature will be to repeat the spectroscopy for smaller fermion masses and various values for 
the diquark source term. 


7. Results at finite density 


At finite density we have no argument that the Pfaffian should remain positive if we add a diquark 
source, but we do no expect a severe sign problem for small baryon densities. Our first simulations 
here were done on an 8 3 x 16 lattice with /3 = 0.96 and K = 0.156. Figure 10 shows the chiral 
condensate, the diquark condensate and the quark number density as a function of small values 
of the chemical potential (below the Silver Blaze transition) and different values for the diquark 
source. Flere we observe that baryon number (and chiral symmetry) is explicitly broken at non¬ 
zero ji such that the quark number density is always non-zero for /./ / 0. We intend to measure the 
diquark spectrum as a function of chemical potential in order to further test our assumptions about 
chiral symmetry breaking at finite density in future. 



P 



P 


Figure 9: Normalized chiral (left) and diquark (right) condensates over the inverse gauge coupling ft. 
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all 


afl 


all 


Figure 10: Chiral condensate (left), diquark condensate (in units of J \, middle) and quark number density 
(right) over chemical potential at T = 0. 


8. Conclusions 

In this contribution to Lattice 2014 we have reviewed our present results for the mass spectrum of 
G2-QCD and its finite density phases at zero temperature. For small quark masses we observe a 
splitting of the spectrum into a pseudo-Goldstone sector, an intermediate bosonic baryon sector, 
and a nucleonic sector. Although the spectrum of CA-QCD is much richer than that of ordinary 
QCD, the results relevant for hadronic physics appear to be quite similar to QCD. In the phase 
diagram at zero temperature one observes several transitions in the baryon or quark number density 
that correlate well with the different mass scales in the baryon spectrum. The transition with 
the smallest value of the chemical potential coincides with the mass of the lightest state carrying 
baryon number, consistent with the Silver Blaze property which is realized in (G-QCD as in two- 
color QCD or QCD at finite isospin chemical potential. The zero-density ground state remains 
unchanged until the quark chemical potential reaches half the mass of the pseudo-Goldstone boson 
which is the scalar diquark as in two-color QCD, a baryon with quark number two. 

We also find evidence that a phase dominated by fermionic baryons exists for quark chemical 
potentials above about 300-600 MeV. In between we found good evidence of a first-order phase 
transition from a phase that is dominated by diquark matter to a phase that is probably dominated 
by ordinary baryonic matter. This transition occurs somewhat below the scale given, in units of the 
quark number, by the nucleon mass as obtained from our vacuum spectroscopy. This might indicate 
a rather large binding energy per nucleon or that, for densities above the Silver Blaze transition, the 
spectrum of the theory depends stronger on the chemical potential than one might expect. 

In order to investigate the chiral properties of the theory in more detail, we have introduced 
diquark sources, in particular, to disentangle the pseudo-Goldstone boson spectrum. This led us 
to simulate Majorana fermions and allowed us to explicitly demonstrate the vacuum realignment 
at zero temperature and chemical potential when the strength of the corrsponding diquark source 
term increases relative to the fermion mass. Moreover, it turned out, that the signal for the finite- 
temperature chiral transition at vanishing net-baryon density is considerably stronger in the diquark 
condensate than it is in the standard chiral condensate. Possible explanations are the explicit break¬ 
ing of chiral symmetry by the Wilson mass that competes with the spontaneous breaking of the 
chiral condensate and the necessity of additive renormalizations. This behaviour is reminiscent 
of simulations with twisted-mass fermions in QCD where the chiral properties of the theory are 
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improved by analogous terms in the fermion action. 
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